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Abstract. We present the most recent searches for lepton-flavor-violating (LFV) T decays in B4B4R. 
We find no evidence of T decaying to three charged leptons or to a charged lepton and a neutral 
meson (A"^, p, ^, K*^K K*^^), and set upper limits on the corresponding branching fractions (BF) 
between 1.8 and 19 x 10"^ at 90% confidence level (CL). 
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THEORETICAL RELEVANCE 

The experimental observation of LFV T and /i decays would provide unambiguous 
evidence of New Physics. In the Standard Model (SM) in fact they proceed through 
diagrams with neutrinos in the loops: the strong GIM suppression, due to the small 
neutrino masses, leads to branching fractions that are below any achievable experimental 
sensitivity. On the other hand, in supersymmetric (SUSY) models, LFV T and /i decays 
can receive significant contributions from diagrams containing SUSY particles in the 
loops, via slepton mixing. Depending on the model and the values of the free parameters 
of the theory (mass spectra and couplings), LFV T decays may have branching fractions 
as high as 10^ (see e.g. [IJ), thus within the reach of the BABAR experiment. 

THE BABAR DATA SAMPLE 

BABAR [2J is a multi-purpose detector operating at the PEP-II e+e collider at SLAC. 
Charged particles' momenta and impact parameters are measured by a tracking system, 
consisting of a silicon strip detector and a gaseous drift chamber in a solenoidal magnetic 
field of 1 .5T. A Tl-doped Csl calorimeter identifies photons and electrons and determines 
their energy. Muons are identified by resistive plate chambers and limited streamer 
tubes installed in the gaps of the iron that contains the solenoid flux return. Kaon/pion 
discrimination is based on the opening angle and photon yield of the Cherenkov light 
emitted in synthetic quartz bars, and on ionization energy loss in the tracking devices. 
Through the years 2000-2008 BABAR has collected almost 1 billion x leptons, produced 
in e+e^^T+T^ events. Here we consider only data taken dX ^ 10.58 GeV (about 
90% of the total T sample). At this center-of-mass (CM) energy the xx cross section, 
~ 0.92 nb, is comparable to that of dimuon production and about 1/5 of the hadronic 
cross section (e^e^ ^qq, q = u,d,s,c,b). The effective cross section for Bhabha events 
where at least one is within the acceptance of the calorimeter is about 50 times larger. 



SELECTION CRITERIA 



We search for the LFV decays T ^ IV^ and IK^, where /(''") = {e,^), V = 

{<l>,p,K*X), (l)^K+K-, p^^Ti+Ti-, K*^^K+7t-, K*^^K-7t+, Kl^Tt+Tt'. We 
look for T-pair events where one (signal) x decays to a fully reconstructed LFV final 
state and the other {tag) T decays to a partially reconstructed SM final state. We consider 
only one-prong tag T decays (BF ^ 85%), except for the eK*^ analysis, where we include 
also three-prong decays (BF ^ 15%). We thus divide the event in two hemispheres in 
the CM frame, by means of a plane perpendicular to the event's thrust axis, and require 
that the 3 charged particles' tracks from the signal T and the tracks (1 or 3) from the tag 
T decay belong to different hemispheres. All tracks must be well reconstructed, within 
the fiducial volume of the detector, and should have zero total charge. Pairs of tracks 
consistent with a photon conversion are removed in order to suppress radiative QED 
backgrounds. The invariant masses of the and daughters are required to be close 
to the mass of the originating mother; the decay vertex should be displaced from the 
interaction point. Tight e/ix/n/K identification criteria, with Ri 1% misidentification 
probability, are applied to the signal T daughters, e/ ix vetoes are applied in some cases 
to the track from 1 -prong tag T decays in order to reduce /i/i and Bhabha backgrounds. 
Finally, since the SM tag T decay contains undetected neutrino(s), missing momentum 
must be different from zero and point inside the detector acceptance. 

From the measured four-momenta of the three decay products of the signal T we deter- 
mine its initial four- momentum and compute the difference Ai? =E'l — ^fsjl between 
the T and the beam energy in the CM, and the difference AM = M^^ — M^^ between 
the energy-constrained invariant mass of the T (after applying the kinematic constraint 
AE — 0) and the T mass world average. In the (AM, A£') plane, signal should peak around 
the origin (with tails at negative AE and positive AM values due to radiation emitted in 
the initial and final state, respectively), while backgrounds should be more uniformly 
distributed. For each decay mode we define, in the (AM, Ai?) plane, a signal region (SR) 
around the origin, and an external sideband (SB) region. After applying tight selection 
requirements in order to reduce the backgrounds as much as possible, we compare the 
number of candidates observed in the SR and the expected background yield, extrapo- 
lated from the number of candidates in the SB. The SR extent and the selection criteria 
are optimized, using large samples of simulated signal and background events in ad- 
dition to background-enriched data control samples, in order to minimize the expected 
branching fraction upper limits (UL) in the background-only hypothesis. 

BACKGROUND EXTRAPOLATION 

Dominant background sources common to all final states are misreconstructed SM xx 

events and random combination of tracks in qq, q = u,d,s events. Important background 
sources are also: Bhabha or dimuon events for the II' I" search and cc events (in particular 
Z)— >y°/v,^^/v decays) for the IV^, IK^ searches. Simulated events and background- 
enriched data control samples (obtained by relaxing or reverting some of the selection 
criteria) are used to extract the shapes of the probability density functions (PDFs) of 
the various background sources in the (AM,AE) plane and their relative abundance. 




FIGURE 1. Distribution of the of the decay tree fit for T — > IK^ candidates selected in data (dots). 
The dashed line and the shaded histogram represent the expected signal and background distributions, 
respectively (arbitrary normalization). 



The normalization of the total background PDF is fitted to the data in the SB, and the 
expected background in the SR is estimated from integration of the PDF in that range. 
After applying the tight selection criteria, only very few background events are expected 
in the SR. Several kinematic distributions and yields in slices of the SB are compared 
between data and simulated events in order to evalutate the reliability of the simulation. 



RESULTS 

For each LFV T decay, the observed number of candidates (A'obs) agrees with the ex- 
pected SM background (A^exp )• Comparing these quantities, taking into account the num- 
ber Nt: of T leptons in the initial sample and the signal selection efficiency e (from 
simulated events), we set upper limits between 1.8 and 19 xlO^^ at 90% CL on the 
corresponding branching fractions. We use either a fully (//'/", IV^ ||6l) or "modified" 
(IKg [|7J) frequentist procedure: in the latter case we compare also the observed and ex- 
pected (in the background-only and background-i-signal hypotheses) distributions of the 
of the geometric and kinematic fit to the whole signal T decay tree. The quoted limits 
include the systematic uncertainties on Nx (~ 1%, from the luminosity uncertainty), e 
(2-9%, dominated by the uncertainty on the particle-identification efficiency) and N^xp 
(typically < 0.3, from the uncertainties in the background PDF shapes and the overall 
background normahzation). The results are summarized in Table[T]and illustrated in Fig- 
ures [T] and |2j In several cases we improve previous experimental bounds. Some regions 
of the parameters' space of some SUSY models (like UJ) are excluded by these results. 
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TABLE 1. Number of T decays (N^), signal efficiency (e), expected background 
candidates (Ngxp) and BF upper limit at 90% CL (BFg^ ), observed candidates (Nobs) 



and branching fraction upper limits (BF^j^j,) for each LFV decay mode. 
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FIGURE 2. (AM, AS) distribution of T ^ ll'l" (top) and T IV° (bottom) candidates selected in data 
(dots). The solid line shows the boundaries of the signal region. The dark and light shading indicates 
contours containing 50% and 90% of the selected MC signal events, respectively. 



